TRANSVERSE POLARIZATION OF QUARK PAIRS 
CREATED IN STRING FRAGMENTATION 



X. ARTRU a , J. CZYZEWSKI 6 

Institut de Physique Nucleaire de Lyon, IN2P3-CNRS et Universite Claude Bernard, 

France 

b Institute of Physics, Jagellonian University, Cracow, Poland 

Abstract 

Classical arguments predict that the quark and the antiquark of a pair created 
during string fragmentation are both transversely polarized in the direction of z x 

, where z is the direction of the pull exerted by the string on the antiquark and 
(— qfjj is the transverse momentum of the quark (antiquark). The existence 
of this effect at the quantum-mechanical level is investigated by considering two 
analogous processes involving the tunnel effect in a strong field : (1) dissociation 
of the positronium atom (2) electron pair creation. In case (1) the positronium 
is taken in the 3 Po state to simulate the vacuum quantum numbers J = ++ . 
Using the nonrelativistic WKB method, the final electron and positron are indeed 
found to be transversely polarized along zxpj_. On the contrary, case (2), treated 
with the Dirac equation, shows no correlation between transverse polarization 
and transverse momentum both when the field is uniform and when it depends on 
z and t. The pair is nevertheless produced in a triplet spin state. The difference 
between these two results and their relevance to transverse spin asymmetry in 
inclusive reactions is discussed. 
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1 Introduction 



The experimental observation of single-spin asymmetries in inclusive hadron produc- 
tion at high energy Q have been tentatively explained by various models : Thomas 
precession 0, P-wave orbitals ||, Regge exchange [|, |5|, semi-classical string mecha- 
nism H [7|]. The asymmetric part of the cross-section is of the form A P± ■ (z x p ± ) 
where z is the collision axis, p ± the transverse momentum of the produced particle, 
P± its transverse polarization (transversity) or that of one of the colliding baryons. 
The "hat" denotes a unitary vector: p ± = p±/\p±\- A similar effect correlating the 
transversity of the leading quark of a jet with the transverse momentum of one of the 
fastest particles was predicted by Collins ||. This effect, if confirmed experimentally, 
would serve as a transverse "quark polarimeter" . It was used in Ref. @ to explain 
single spin asymmetry in inclusive meson production. 



In this paper we start from a popular string hadronization picture |5], |IU], [TJJ, in 



which quark-antiquark pairs are produced from the string by a tunneling mechanism 
analogous to the Schwinger mechanism for pair creation in a strong homogeneous elec- 
tric field [12-20]. This picture accounts rather well for the exponential cut-off in p ± , the 
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Figure 1: Semi-classical string mechanism of quark polarization. The orbital angular 
momentum of the qq pair is compensated by the spin of q and q, thereby causing the 
correlation between spin and transverse momentum of the quark and the antiquark. 



relative suppression of strange quarks and the almost complete suppression of heavy 
quarks. Let us recall the semi-classical arguments |J for a transverse polarization of 
the quark and the antiquark (see Fig. 1) : 

- the quark and the antiquark come from a pair fluctuation like those which occur in 
ordinary vacuum. At zero separation, the pair has zero total energy-momentum. 
In particular, quark and antiquark have opposite transverse momentum q ± and 
q ± . In the vacuum case, the pair stays virtual and disappears after a time of 
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the order of the quark Compton wave length. In the string case, the linear mass 
density k ~ 1 GeV/fm of the string is converted into energy of the pair, which 
becomes real at a longitudinal separation d = 2Ej_/k, where E± = [m 2 + q\) 1 ^ 2 
is the quark (or antiquark) transverse energy. 

- The orbital angular momentum 

L = d z x q ± (1) 

is compensated, at least partly, by the spins of the quark and the antiquark. 
Assuming equal polarization P = 2 < s q > = 2 < Sq > for the quark and the 
antiquark, we have therefore 

P = -Lf(L), (2) 

where f(L) is a reduction factor insuring that \P\ is smaller than unity, e.g., 
f(L) = (1 + L)-\ 

To summarize, the polarizations of the quark and of the antiquark are of the form 

P q = P q = -A(q ± ) zxq ± , (3) 

where 

A(q ± ) = Lf(L) <min{L,l} (4) 
is the analysing power of the mechanism and 

L = 2K- 1 q ± (m 2 + q 2 ± ) 1/2 . (5) 

The compensation between L and s q + s q is further motivated by the phenomenological 
success of the u3 P " model of quark pair creation in hadronic decay [pi]. This model 



assumes that the (q q) pair is created in the 3 P spin state, therefore having the vacuum 
quantum numbers J PC = ++ (by contrast, a model in which the pair comes from one 
gluon gives 

The string, as well as the constant electric field, is not rotationally invariant and 
therefore the total angular momentum J = L + s q + Sq is not conserved during tun- 
neling. This is a weak point of the classical model reviewed above. C- and P- quan- 
tum numbers are also not separately conserved. Therefore it is desirable to check the 
existence of the transverse polarization effect in a true quantum-mechanical model. 
Although too difficult at present in the string theory, this is quite possible for the 
problem of pair creation in strong homogeneous field. To begin with (in section 2) 
we will consider the nonrelativistic process of the dissociation of a positronium atom, 
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which we assume to be in the 3 Pq state. The relativistic case of pair creation in a field 
S(t, z) parallel to the z axis (independent of x and y but not necessarily on t and z) will 
be considered in section 3, using the Dirac hole theory for the positron. The different 
results obtained in these two cases will be discussed in Section 4. 



2 Positronium dissociation 

The nonrelativistic e + e~ system in a constant electric field is governed by the Hamil- 
tonian 

Oi pt 

1 + e 5 ■ (r_ — r H 
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+ K r + V c -Fz. (6) 



The motion of the barycentre can be separated from the relative motion and from 
now on we will consider only the latter, governed by the three last terms of Eq. (^|). 
r = r + — r_ and p = (p , — p_)/2 are the relative position and momentum. We take 
F = eS z > 0. The potential V c — Fz is shown in Fig. 2. Note that the Hamiltonian 
is spin-independent, therefore cannot produce spin effects by itself. However, we will 
assume that the initial state of the pair is a 3 Pq positronium (corresponding to the 
vacuum quantum numbers) : 

$(r) = f(r) [ Yr \r ) | + 1) - Y?(r) |0> + n +1 (f ) | - 1) ] , (7) 

where the kets denote the three different triplet spin states, 

| + i) = ITT), |0) = (|TI) + |IT))/v / 2, |-i) = |||). (8) 

In this way, the orbital motion and the spin are entangled. This state is a bound 
eigenstate of K r + V c with energy —B. After turning on the external electric field, the 
relative wave function will eventually migrate toward z = +oo by tunnel effect, which 
means that the positron runs toward z = +oo and the electron runs toward z = — oo. 
The pair remains in the spin-triplet subspace. We chose x as the spin quantization 
axis and are interested in the relative probabilities to obtain the different final spin 
states \S X ), with S x = s+ + s~ , for a given final transverse momentum q ± = qy. The 
corresponding asymptotic state is an eigenstate of K r — Fz, again with energy —B, 
and its wave function is 

*(r) = V(r) 1^) = g(z - z t ) \S X ) (9) 
where g is the Airy function, solution of the one-dimensional Schrodinger equation 

'p 2 Jm-F{z-z t )} g(z) = (10) 
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Figure 2: Positronium dissociation in constant electric field. The top curve is the 
superposition of the Coulomb potential and the external electric potential. The clas- 
sically allowed region for a given transverse momentum q is limited by the horizontal 
dashed line. The bottom curves are the p wave function of the positronium, and 
the wave function of the free solution, 'J (restricted to the z axis) . The overlap of <I> 
and ^ is responsible for the tunneling. 



and zt = (B + q 2 /m)/F is the classical turning point (see Fig. 2). Here we give a 
heuristic proof and an estimation^ of the spin asymmetry : 

- We assume that the tunneling length z t ~ ma 2 /(16F) is much larger than the 
radius ~ 8/ (ma) of the bound state. Near the bound state, we can use the WKB 
approximation for g : 

g(z - zt) ~ e A <*-*> (11) 

with A = (mB + q 2 ) 1 / 2 . 

- Near the origin, ip(r) can be expanded in partial waves : 

il>(r) ~ (&' r = AnJ^i 1 3i(pr) Y f m {P) Y i m {r) (12) 

l,m 

1 A more rigourous treatment could be done using the method of Landau & Lifshitz (Quantum 
Mechanics p2j), for hydrogen dissociation in a strong electric field. 
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with p = qy — i\z, p = i(\ 2 — q 2 ) 1 ^ 2 , p = p/p. We assume that tunneling couples 
mainly the components of $ and \I/ with the same YJ m (f), and the tunneling 
amplitude is proportional to the coefficient of this harmonic (this is intuitive 
if we consider the inverse process of trapping an initially free particle into the 
Coulomb potential well). The I = 1 terms of ip are proportional to 



j x {pr) \p 2 \ 1/2 [ -{p y - ip z ) Y+\r) + (p y + ip z ) Y 1 l \ 



r 



= hipr) (A 2 - q 2 r 1/2 [ (A - g) Y?\*) + (A + q) Y^r) ] , (13) 

Comparing with Eq. (0) we find that the tunneling amplitudes squared are in the ratio 

|T(5, = +1)| 2 : |T(5, = 0)| 2 : \T(S X = -1)\ 2 = \X + q\ 2 : : |A - q\ 2 (14) 

Note the vanishing of T(S X = 0). It happens because the second term of $ is odd in 
x and cannot tunnel to which is even in x (orbital x-parity is a symmetry of the 
problem). The polarization of the electron and the positron are equal and given by 

\T( + 1)\ 2 - \T(-1)\ 2 2 y /mBT~q 2 
^ ~ ir(+l)| 2 + 35 _ mB + 2q 2 * q± 

We see that the polarization of the created particle is of the form @, (||) and has the 
same sign as predicted by the classical string arguments. Classical trajectories leading 
to the positronium dissociation shown in Fig. 3 explain this fact intuitively. 




fa voured case unfa voured case 

Figure 3: Classical trajectories of the positronium dissociation for the two cases 
L x — +1 and L x = — 1. The dashed lines are (classically forbidden) tunneling trajec- 
tories. 



3 Pair creation in strong field 

In a static constant electric field, electron-positron pairs are created spontaneously in 
vacuum, the positron running in the field direction and the electron in the opposite 
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direction. In the Dirac hole theory, this process is interpreted as tunneling of the 
electron from the Dirac sea in one half-space to the upper continuum in the opposite 



half-space, without changing its total energy, as shown in Fig. 4 [12]. We will study 



Zl z 



V+m 
V 

V-m 

Figure 4: The Dirac sea distorted by constant electric field between and Z\. A 
negative-energy electron on the left-hand side can reach the upper continuum of the 
right-hand side by tunneling trough the forbidden band, becoming physical electron. 
The hole created on the left-hand side is the physical positron. The dashed line 
represent the energy of the tunneling wave function. 

the more general case of a time- and z-dependent field S(t, z) parallel to z and consider 
Dirac wave functions of definite transverse momentum parallel to the y axis : p ± = 
(Px,P y ) = (0,g). Discarding the trivial y-dependence in exp(iqy), the Dirac equation 
reads 

[id t + eA - a z (id z + eA z ) - qa y - m(3] tp(t, z) = . (16) 
Calculations will be simpler using the light-cone coordinates 

d={t-z)/2 (17) 

dt = d t -d z (18) 
As = A - A 3 = A + A 3 (19) 
Furthermore we choose the spinorial representations of the Dirac matrices 
/<r* \ /OA , / a i \ 

where a 1 are the Pauli matrices. The Dirac equation becomes 

[a v (id v + eA v ) + (id^ + eA^) - qa y - m(3} ip(r],£) = 




V = (t + z)/2, 

d v = d t + d z , 
A V = A + A 3 = A°-A 3 , 
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V 



^,0 = o 



(21) 



/ id v + eA v iq — m 

— iq id^ + eA% — m 

—m id^ + — iq 

— m iq id n + eA^ / 

This equation is invariant under the transformation ^i « t/>4, ^ ^ ^3, which is 
performed by the matrix (3H X = ^j 5 . This matrix commutes with the Hamiltonian 
and has eigenvalues ±1. For a particle at rest, (3H X — 2s x . For a particle with 
nonzero p y and p z , it is the x-component of the transversity operator. For the states 
under consideration, the (3 Y^ x transformation is equivalent to the parity about the (y,z) 
plane, 



P 



-iirJx 



i/z 



P 



-i S x X P X intrinsic x P, 



orbital 



(22) 



Since p x = 0, P™" = 1 and (3 Y> x is equivalent to P yz (up to a phase factor, depending 
on the choice of the intrinsic parity) . The (3 S x invariance comes therefore from the 
symmetry of the problem about the (y,z) plane. 



Taking the transversity eigenstates 



G 
G 



^ = 7=2 



( F \ 

G 

-G 
\-F) 



(23) 



we come to coupled differential equations 



(id n + eA v ) F = (±m - iq) G (24) 

(z<% + eAt.) G = (±m + iq) F (25) 

where ± is the sign of the transversity. From these we get the second order differential 
equations 

[ (idt + eA^) (id v + eA v ) - M 2 ] F = (26) 

[ (id v + eA v ) (id 6 + eA^ - M 2 } G = , (27) 

where M = (m 2 + g 2 ) 1 / 2 is the "transverse energy" of the electron. Note that Eqs. (p6|) 
and ( p7|) depend only on the transverse energy, not on the transversity. We can infer 
that there is no correlation between transverse spin and transverse momentum, con- 
trarily to the case of positronium dissociation. Let us check this result more carefully. 
Denoting by F mA and G m ^ the solution of Eqs. (|2l|]2"5| ) and setting m + iq = Me™, 



we have 



m,q 



F M ,o e 



=F«o/2 



n — n C ±W 2 



(2* 
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These equations tell how the solution transforms under a "rotation" in the (m,q) plane 
(which leaves M invariant). The current 4- vector (J°, J 1 ) = (ip'tpjip'ofip) is given by 

jo = |F| 2 + |Gf, J x = 0, J Z = \F\ 2 -\G\ 2 (29) 
jy= 2lm(F*G) = ^R e [F*(td r] + eA v )F]±^d v \F\ 2 . (30) 



The last form of J y has been obtained using Eq. (|24|) . The first part is the "convection" 
term, the second part is the "magnetization" term. J°, J z and the convection part 
of J y do not depend on the transversity. This confirms the absence of transverse spin 
effect in the Schwinger mechanism of pair creation. The magnetization term of J y 
depends on transversity but is located at the edges of the wave packet (we can replace 
9,j|F| 2 by (do J z + d z J°)/2, using current conservation) and is not observable by a 
macroscopic detector. 

To fix the idea, let us consider a homogeneous field £ = £ z confined in the region 
< £ < £i. This field corresponds to a capacitor moving with light velocity. We use 
the null-plane gauge 

( if £ < 
A^ = 0, A v =<2££ if0<£<£i (31) 

[2^£ if a < e 

At fixed light-cone momentum p v = p +p 3 = p° —p z , we have the following solutions : 

G = e- ip ^~ iAp , F = ±m ~ tq G, (£ < 0) (32) 

Pv 

G = e-*» f^ + lV" , F=^-^-G, (0<£<£i) (33) 
V Pv J Pv + 2 K £ 



G= U +1 ) exp 



M 2 - 

■ T.jri £ ~~ £1 

tp v r] — iM 



±m -' q G, ({,<£) 



p v + 2 k£i 



p v + 2 k£i 

(34) 

k = e£ is the electric force and P v = p v and P^ = p v + 2 k£i are the initial and final 
"mechanical" (gauge invariant) light-cone momenta (P M = p M + eA M ). Electrons from 
the Dirac sea at £ < having light-cone momentum P v in the range [— 2/c£i, 0] become 
physical electrons (P' > 0) at 

£ c = -P„/ (2/c) (35) 

The electron flux going through the hyperplane £ = constant being proportional to 
J ? = (J° — J z )/2 = \G\ 2 , the tunneling probability is 

^(^£>£ c )/^(^e<£e)=e- M2/K . (36) 
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This result is clearly independent on spin. 

Some remarks have to be made concerning the above calculation: 

- The last result is obtained giving a small positive imaginary part to P v . This 
corresponds to the physical condition that the field does not interact with the 
wave at t — — oo. 

- In such a field, the created electron escape the field region (at £ = £1) but not 
the positron. This can be seen from their classical trajectories in the field region 
shown in Fig. 5 : 

rj = Vo ± M/ {2k) e r , £ = £ c =F M/ (2k) e~ r , y = y + qr/K (37) 

where 770 and yo are f ree parameters and r = ± (k/M) x (proper time) is the 
rapidity of the e ± . 




Figure 5: Classical trajectories of the electron and the positron described by Eq. (37) 
created in the field Eq. (31) confined in the region bounded by the two solid diagonal 
lines. The electron escapes from the field region while the positron remains in the 
field forever. 

- J v — \F\ 2 becomes infinite at £ = £ c . The current looks like a "jet stream". It is 
due to the deflection of the incoming flux by the field during infinite time. 
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4 Discussion 



After obtaining the positive result with the positronium model, the absence of trans- 
verse polarization in the Schwinger mechanism was rather unexpected. This absence 
does not happen due to standard discrete symmetries like C, P and T but due to 
invariance with respect to the particular transformation (f2§|). 

In spite of the difference between the positronium dissociation and the Schwinger 
mechanism, both models predict that the electron and the positron have equal spin 
components s+ and s~ along the x axis : 

sis- = +1/4 (38) 

This property is built-in in the 3 Pq positronium model. For the Schwinger mechanism, 
both s+ and s~ are equal to ^/3T> X . Indeed, a positron at rest has s+ = — since 
it is a "hole", and (3 = — 1 because the corresponding (unoccupied) state has nega- 
tive energy. Eq. (^) imply that the pair is in a triplet state also for the Schwinger 
mechanism. 

The difference between the two models may be connected with their different chiral 
properties, which appear in the m —>■ limit : 3 P positronium dissociation is similar 
to the decay of a ++ particle into two fermions, in which the fermions necessarily have 
opposite chirality (equal helicity). On the contrary, the Schwinger mechanism involves 
only vector interactions, therefore e + and e~ must have the same chirality (opposite 
helicity). The difference is particularly important at zero transverse momentum, where 
conservation of angular momentum J z requires s~ = —s^. Positronium dissociation 
is allowed for m = q± = 0, whereas pair creation (with back-to-back e + and e~) is 
forbidden?]. 

The absence of correlation between spin and transverse momentum in the Schwinger 
mechanism does not preclude such correlations for qq pairs created during string break- 
ing. There are many effects of string breaking not included in the Schwinger mech- 
anism. First of all, the chromoelectric field between the quark and the antiquark is 
totally screened after their creation, unlike in the Schwinger process where the field 
extends everywhere all the time^ . Secondly, the field of the QCD string is confined to 

2 Eq.([36|) seems to allow pair creation at m = q± = 0. However, this is a too special case where 
the field (|3l| ) occupies an infinite domain of the (z,t) plane. In fact, for m = q± = 0, the functions 
F and G decouple (see Eqs. ^6 27) and the left- and right-moving currents = \G\ 2 and J n = \F\ 2 
are separately conserved. If the field domain is finite in the (z,t) plane, e + e _ pairs are produced at 
m = q± = 0, but with e + and e~ going in the same direction. e + and e~ going in opposite direction 
belong to independant pairs. 

3 Let us mention however Ref. [19| where the screening effect has been considered for scalar QED. 
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a thin tube. One way to simulate this fact in QED is to impose the MIT-bag boundary 
conditions for the electron in the transverse coordinates [IT], In our problem, we 
cannot use this method because we need a well-defined transverse momentum. Impor- 
tant effects may also come from the "transverse inertia" of the string, because part of 
it must follow the transverse motion of the quark. Finally, the Schwinger mechanism 
does not include the final state interactions which recombine the quarks from different 
pairs to form hadrons and resonances. Resonances probably play a major role in single 
spin asymmetries || |2^, because the latter come from the interference between 
different spin amplitudes having different phases. 

To summarize, the interplay of spin and transverse momentum in pair creation is 
a subtle phenomenon and we cannot conclude from the simple model presented above 
whether the correlation exists or not. 



Acknowledgements 



We acknowledge the financial support from the IN2P3-Poland scientific exchange pro- 
gramme (collaboration no. 91-62). J.C. has been also supported by the Polish State 
Committee for Scientific Research (KBN) grant no. 2 P03B 086 14 and by the Polish- 
German Collaboration Foundation grant FWPN no. 1441/LN/94 during complition of 
this work. 



11 



References 

[1] E704 Coll., A. Bravar et al, Phys. Rev. Lett. 77, 2626 (1996); E704 Coll., 
D.L. Adams et al, Z. Phys. C56, 181 (1992); E704 Coll., D.L. Adams et al, 
Phys. Lett. B264, 462 (1991); E704 Coll., D.L. Adams et al, Phys. Lett. B261, 
201 (1991); V.D. Apokin et al, Phys. Lett. B243, 461 (1990); S. Saroff et al, 
Phys. Rev. Lett. 64, 995 (1990); M.S. Amaglobeli et al, Sov. J. Nucl. Phys. 50, 
432 (1989); E704 Coll., B.E. Bonner et al, Phys. Rev. Lett. 61, 1918 (1988); 
J. Antille et al, Phys. Lett. B94, 523 (1980). 

[2] T. DeGrand, H.I. Miettinen, Phys. Rev. D24, 2419 (1981); T. DeGrand, J. 
Markkanen, H.I. Miettinen, Phys. Rev. D32, 2445 (1985). 

[3] C. Boros, Liang Zuo-tang, Phys. Rev. D53, R2279 (1996). 

[4] J. Soffer, N.A. Tornqvist, Phys. Rev. Lett. 68, 907 (1992). 

[5] C. Barni, G. Preparata, PG. Ratcliffe, Phys. Lett. B296 251 (1992). 

[6] B. Andersson et al, Phys. Lett. B85, 417 (1979); B. Andersson, G. Gustafson, 
G. Ingelman, T. Sjostrand, Phys. Rep. 97 31 (1983). 

[7] X. Artru, J. Czyzewski, H. Yabuki, Z. Phys. C73, 527 (1997). 

[8] J. Collins, Nucl. Phys. B396, 161 (1993); J. Collins, S.F. Heppelmann and 
G.A. Ladinsky, Nucl. Phys. B420, 565 (1994). 

[9] A. Casher, H. Neuberger, S. Nussinov, Phys. Rev. D20, 179 (1979). 

[10] E.G. Gurvich, Phys. Lett. B87, 386 (1979). 

[11] N.K. Glendenning, T. Matsui, Phys. Rev. D28, 2890 (1983). 

[12] J. Schwinger, Phys. Rev. 82, 664 (1951); 

[13] E. Brezin, C. Itzykson, Phys. Rev. D2, 1191 (1970). 

[14] R.C. Wang, C.Y. Wong, Phys. Rev. D38, 348 (1988); C.Y. Wong, R.C. Wang, 
J.S. Wu, Phys. Rev. D51, 3940 (1995). 

[15] C. Martin, D. Vautherin, Phys. Rev. D38, 3593 (1988). 

[16] M. Herrmann, J. Knoll, W. Greiner, Phys. Lett. B234, 437 (1990). 



12 



[17] Th. Schonfeld, A. Schafer, B. Muller, K. Sailer, J. Reinhardt, W. Greiner, Phys. 
Lett. B247, 5 (1990); K. Sailer, Th. Schdnfeld, A. Schafer, B. Miiller, W. Greiner, 
Phys. Lett. B240, 381 (1990); K. Sailer, Z. Hornyaak, A. Schafer, W. Greiner, 
Phys. Lett. B287, 349 (1992). 

[18] H.-P. Pavel, D.M. Brink, Z. Phys. C51, 119 (1991). 

[19] Y. Kluger, J.M. Eisenberg, B. Svetitsky, F. Cooper, E. Mottola, Phys. Rev. Lett. 
67, 2427 (1991). 

[20] S.A. Smolyanskii, G. Ropke, S. Schmidt, D. Blaschke, V.D. Toneev, A.V. Pro- 
zorkevich, Preprint GSI 97-72 (Dec. 1997), |hep-ph/ 971 2377 . 

[21] L. Micu, Nucl. Phys. BIO, 521 (1969); A. Le Yaouanc, L. Oliver, O. Pene, J.- 
C. Raynal, Phys. Rev. D8, 2223 (1973), Hadron Transition in the Quark Model 
(Gordon Breach, London, 1988). 

[22] L. Landau, E. Lifshitz, Course of Theoretical Physics, Vol. 3, Quantum Mechanics 
(third ed., 1975). 

[23] J. Collins, G.A. Ladinsky, PSU/TH/114 (Dec. 1994); |hep-ph/9411444 . 
[24] R.L. Jaffe, Xuemin Jin, Jian Tang, Phys. Rev. Lett. 80, 1166 (1998). 



13 



